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Abstract 

A bound on the error introduced by truncating a quantum addition 
is given. This bound shows that only a few controhed rotation gates 
win be necessary to get a rehable computation. 



Introduction 

Draper introduced an addition algorithm for a quantum computer that 
uses only Hadamard and controlled rotation gates. For an n-bit addition, 
his algorithm requires only 2n qubits, O(n^) gates in its full implementation, 
and allows for parallelization. He also cited Barenco, et. al., to appeal that 
an approximate quantum Fourier transform may be more accurate than a 
full quantum Fourier transform ||2|. In this paper, we investigate the error 
introduced by truncating the algorithm by removing all gates with tolerance 
less than a certain threshold, assuming that no decoherence is introduced into 
the computation. This work shows that the probability of error decreases 
rapidly as one adds more gates, so that reliable computation can be achieved 
even with relatively few controlled rotation gates used in the approximate 
algorithm. The reader is assumed to be familiar with Draper's paper, as this 
work will draw heavily from it as a source of ideas. 

From a high level, the full machinery of Draper's algorithm goes per the 
following diagram, where \a) and \b) are both n-bit quantum registers, and 
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QFTgn is the quantum Fourier transform on n qubits. 

\b) ^ \b) \b) — \b) 

\a) — > QFTan \a) — > QFTgn \a + b) — > \a + b) ^ ^ 

Here a + b represents the n-long bit string formed by the sum of a and b 
considered as numbers modulo 2"-, not the vector sum of \a) and \b). 

The QFT algorithm used in this computation is a refinement due to Cop- 
persmith that shows how the QFT can be implemented using only Hadamard 
and controlled rotation gates [0. A controlled rotation gate is given by the 
function Rg, which is the action 

Re : \a) ® \b) ^ e'''''' \a) ® \b) , 

where a and b are either or 1. We are supposing that we cannot perform a 
rotation Re for an angle \9\ < 2n/2'^ for some k because such a gate may cost 
too much to manufacture, may not be reliable, or the computation may be 
faster because of a fewer number of gates. The important number here is k 
because we are calculating the error introduced by truncating a computation 
by eliminating all gates Re, where \6\ < 27r/2'^. We will call the maximal k 
the threshold for our computation. We would like to replace Diagram || 
with the following diagram, where A is the approximate quantum Fourier 
transform on n qubits with k as the threshold, and represents controlled 
rotations conditioned on b. 

\b) \b) \b) \b) 

\a) — > A\a) — > ^6v4|a) — > A'^'^bAla) 

We want to ensure that the quantity 

\{a + b\A^<i/bA\a)\\ 

which is the probability that the the approximate addition algorithm yields 
the same result as the full implementation, is close to unity for a given 
threshold k. 

Notation and Background 

A review of notation is needed at this point. We will follow Draper's notation 
by defining e{t) = e^'^**, so that e{x + y) = e{x)e{y), and e{x) = 1 if and only 
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if X is an integer. We also define 

U.a„a„_i . . .a^ = — h • • ■ H , 

where each is either or 1. Note that this is equivalent to a/2", where 

a = a^ttn-i . . . Oi = a„2"~^ + a„_i2"~^ + ■ ■ ■ + ai2^. 

Then if a is an ra-bit integer with n > k, we have 

e(a/2^') = e(a„2"-^-i + a„_i2"-^-2 + ■ ■ • + a^+i + ak/2 + ■■■ + 
= e{ak/2 + ak-i/2^ + ■■■ + ai/2'') 
= e(0.afcafc_i . . . ai). 

The quantum Fourier transform on an iV- dimensional space evolves the 
basis state vector \a) to the state 



7V-1 

QFT^ : l«) - ^ E 1^') ' (2) 



N-l 

N . „ 

where the orthonormal basis vectors of the space are enumerated |0), |1), . . . , 
— 1). If is composite, say = MK for some integers M and K, then 
we can break up the quantum Fourier transform into the tensor product of 
at least two distinct unentangled parts. By the division algorithm, we can 
let j = mK + k, where < m < M, and < k < K. We assume that we are 
using the positional notation with radix K, so that \j) = |m) (S) \k). After 
substituting for j and A^, the sum in Equation || becomes 

M~1K~1 

QFT^ \a) = -jL= Yl e2--(™^+^)/^^ \m) ® \k) 

^ m=0 k=0 

AI-1 K-1 



(QFT,,|a))®(-^ge^™^/^|A:)). 



(3) 
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Equation |^ deserves some commentary, in that we seem to be confusing math- 
ematical notation for actual physical reality. The content of this equation is 
that we can take the A^-dimensional space An and decompose it as a tensor 
product Am ® Ak- We then get the QFT of \a mod M) in Am, and tensor 
that with the result of the operation specified in Ak, with the final result of 
the computation being the QFT of \a) in Ajy. Hopefully this sloppiness will 
not cause to much confusion to the reader, as it is a mathematical artifice 
too show that the quantum Fourier transform is not entangled. Now we as- 
sume that = 2" is a power of 2. We let K = 2, and M = 2"~^, and then 
substitute M and K into the Equation ^ above, finding the equation 

QFT2„ \a) = (QFT2,.-i |a)) ® ^(|0) + e(a/2'^) |1)). (4) 

V2 

We continue inductively, eventually arriving at: 

QFT^n \a) = ^(|0) + e(a/2) |1)) ® (|0) + eia/2') |1)) ® ■ • ■ 

v2" [0) 

® (|0)+e(a/2")|l)). 

This formula allows an implementation of the n-bit quantum Fourier trans- 
form using a total of n Hadamard gates, n(n — 1)/2 controlled rotation gates, 
and [n/2\ swap gates, as shown in Coppersmith's paper [^. 

The Approximate QFT 

The approximate quantum Fourier transform (which I call A throughout 
this paper) is implemented by deleting all controlled rotation gates with 
a tolerance finer than a certain threshold. If the threshold is k, then the 
approximate QFT on n qubits is given by 

A \a) = ^(|0) + e(a/2) |1)) ® (|0) + e{a/2') |1)) 

® ■ ■ ■ ® (|0) + e(a/2^) |1)) ® (|0) + e{0.ak+iak ... 02) |1)) (6) 
O (|0) +e(0.afc+2afc...a3) |1)) O--- 
(g) (|0) + e(0.a„a„_i . . .a„_fc+i) 

Written this way, the approximate QFT can be put into a useful form for 
analysis. First, we need to introduce some new notation. Define Tr^(a), 
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representing truncation, as being a function that takes the number 

a = a„2"-^ + a„_i2"-2 + . . . + ai2° 

and nonnegative integers k and m, where m < n, and yields the number 

Tr^(a) = 0. OO^a^a^.i . . . ai, 
A; O's 

so that this is the m least significant digits of a multiplied by 2^^™+'^). Note 
that 

< Tr^(a) + Tr^b) - Tr^(a + 6) < 2-^ (7) 

which captures the fact that carry bits propogate through the addition. The 
error is bounded above because the only possible non-zero difference between 
the sum of the truncations and the truncation of the sum is the number 

0. pOOOOO...Q l = 2-^ 

A; - 1 O's 

From (0), it is easy to make a inductive argument to prove that 



< (ETr-(a,)) - Trr (X^a,) < (n - 1) ■ 2"^ 

Also, it is easy to see that 

e{a/2"^) = 6(0.0^0^-1 ... ai) 

= e{0.amam-i ■ ■ ■ am-k+i + Tr^"''(a)). 

We can therefore write the approximate quantum Fourier transform as 
A \a) = 4^(10) + e(a/2) |1)) ® (|0) + e(a/2^) |1)) 

V 2" 



(8) 



+ e{a/2^) |1)) ® (|0) + e(a/2^+i - Tr^(a)) |1)) (g) 
+ e(a/2^+2-Tr^(a))|l))®--- 
+ e(a/2"-Trr^(a))|l)). 
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Rotation by b 



We have seen how to compute the approximate quantum Fourier transform 
of a state vector \a). At this point, we want to see the action of controlled 
rotations conditioned on bits from b. If we take the last qubit in Equation 
we see that after the application of the approximate QFT, its unnormalized 
state is 

\4>n) = |0) + e(0.ana„_i . . . an-k+i) |1) 
and is rotated to the state 

\(pn) = |0) + e(0.ana„_i . . . a„_A:+i + 0.6„6„_i . . . bn-k+i) |1) 
by the action of \l'f, This state is equivalent to 

^6 |0n) = |0) + e((a + 6)/2" - TT^^a) - Trr^(6)) |1) . 

Note that the correct state, i.e., if we were to perform an approximate QFT 
on the vector \a + b), is 

l^n) = |0) + e((a + 6)/2" - Trr'=(a + b)) |1) . 

Similar statements hold for the other qubits. We discover the error intro- 
duced by using the approximate transform by rewriting this state as 

^6 |0n) = m + {e((a + 6)/2") [e(- Trr'^(a) - Trr'=(6)) 
-e{-T4-\a + b))]\l)]. 

Then \ipn) will be evolved back to the correct answer of a+b by the application 
of the inverse approximate QFT, . The remainder of the state vector 
enclosed in braces is an error term. There is also a normalization factor of 
1 / V2 that I have not included. 

After the application of and the approximate QFT, the a register has 
been evolved to a state 

^bA \a) = l^i) ® IV'2) ® ■ ■ ■ ® l^fe) ® il^Pk+i) + \lk+i)) .^Q. 

® + |7fe+2)) ® ■ ■ • ® (l^n) + |7n)), 
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where |7j) represents the error on qubit j. If we expand this expression 
out, we have 2"'"'^ terms. There will be f""'^) terms with m errors for each 
< m < n — k. We have the identity that 



ze** dt. 



So 



< 



ie** dt 

I 

\ie''\dt 

e-a\ 



(11) 



for any 6 and a. Combining this with the inequality (|^) yields that, for every 
j such that n — k < j < n, 



|7,) \\<nV2-2-'. 



We also know that 



\a) ® \b) 



l«) II ■ II \b) 



and that || lipi) || = 1 for every 1 < i < n; therefore, the effect of having 
m errors in a term is that the magnitude of the error on that term is equal 
to the product of the magnitudes of each of the m errors. Since A is unitary, 
the magnitude of the error vector is preserved when we apply Using 
the triangle inequality (|| X]i^«ll — J2i 11^*11) to bound the error in the state 



given in Equation |T0|, we see that the magnitude of the error must be bounded 
above by 



m=l ^ ^ 



\—k\n—k 



The probability of error is this number squared. When we do a first-order 
expansion of this expression, we get 



TT{n -k)-2-^ -^2 



(12) 
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as a rough estimate for the bound. Thus k does not have to be very large 
compared to n to maintain a very small probability of error in the compu- 



tation. For example, if n = 1000, and k = 30, then the bound given in (|12|) 
is approximately 4.0 x 10~^; the probability of getting an error is smaller 
than 1.6x10"^^. One can use (H) to prove that if we perform m /c-truncated 
additions — i.e., add m + 1 numbers — the magnitude of the error vector is 
bounded above by 

{1 + 7rV2 ■ 2-'' ■ my-'' - 1, 

which also goes quickly to zero, and thus implies that there is a high prob- 
ability of getting the correct answer when one uses a sufficient number of 
controlled rotation gates. 



Summary 

We have derived a bound on the error introduced by truncating in Draper's 
quantum addition algorithm. This bound goes quickly toward zero, and 
therefore promises that only a moderate number of rotation gates are required 
to ensure reliable addition if we use Draper's algorithm in connection with 
reliable quantum computational gates. 



References 

[1] T. Draper. "Addition on a Quantum Computer". |quant-ph/000803^ . 

[2] A. Barenco, A. Ekert, K. Suominen, P. Torma. "Approximate Quantum 
Fourier Transform and Decoherence" . [quant -ph/ 960 101^ 

[3] D. Coppersmith. "An Approximate Fourier Transform Useful in Quantum 
Factoring". Technical Report, IBM, 1994. IBM Research Report 19642. 



8 



